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We present a study of the spectral properties of photon pairs generated through the process of 
spontaneous four wave mixing (SFWM) in single mode fiber. Our analysis assumes narrowband 
pumps, which are allowed to be frequency-degenerate or non-degenerate. Based on this analysis, we 
derive conditions on the pump frequencies and on the fiber dispersion parameters which guarantee 
the generation of ultra-broadband photon pairs. Such photon pairs are characterized by: i) a very 
large degree of entanglement, and ii) a very high degree of temporal synchronization between the 
signal and idler photons. Through a numerical exercise, we find that the use of photonic crystal 
fiber (PCF) facilitates the fulfilment of the conditions for ultra-broadband photon pair generation; 
in particular, the spectral region in which emission occurs can be adjusted to particular needs 
through an appropriate choice of the PCF parameters. In addition, we present a novel quantum 
interference effect, resulting from indistinguishable pathways to the same outcome, which can occur 
when pumping a SFWM source with multiple spectral lines. 

PACS numbers: 42.50.-p, 03.65.Ud, 42.65.Hw 



I. INTRODUCTION 

Quantum-enhanced technologies require photon pairs 
with specific properties. Thus, for example, factorabil- 
ity is required for the generation of pure heralded sin- 
gle photons, a crucial resource for linear optical quan- 
tum computation [T]. At the opposite extreme, spon- 
taneous parametric processes permit the generation of 
photon pairs with a remarkably large degree of entangle- 
ment. Such photon pairs are characterized by a large 
Schmidt number [2], which implies that emission takes 
place into a large number of independent frequency or 
transverse momentum signal and idler mode pairs. A 
state with these characteristics leads to a large mutual 
information, which quantifies the information which two 
parties can in principle share by virtue of the entangle- 
ment present [51. Two-photon states with a large degree 
of entanglement may lead to exciting applications, such 
as large-alphabet quantum key distribution |3], quantum- 
enhanced two photon absorption [5], and teleportation of 
single-photon wavepackets[5]. 

Most research on photon pair generation has relied on 
the process of spontaneous parametric downconversion 
(PDC) in second-order non-linear crystals[7j. Recently, 
the process of spontaneous four-wave-mixing (SFWM) 
in optical fiber, relying on a third-order non-linearity, 
has emerged as a useful alternative, with some clear ad- 
vantages. Indeed, fibers permit an essentially unlimited 
interaction length; in addition, the generated flux is pro- 
portional to the square of the pump power, instead of 
the linear pump-power dependence observed for PDC. 
This leads to the possibility of remarkably bright sources. 
While some of the first SFWM experiments were carried 
out with standard telecom fiber [H], recent experiments 



have exploited the greater flexibility conferred by pho- 
tonic crystal flbers (PCF) [9]. Thus, for example, the zero 
group velocity dispersion frequencies may be selected by 
careful choice of the flber design parameters, leading to 
the ability to also select the generation frequencies. In 
previous work, we have shown [TO] that SFWM in PCF's 
leads to the ability to engineer the spectral entanglement 
properties of photon pairs. 

We are particulary interested in spontaneous paramet- 
ric processes constrained to a single transverse mode, e.g. 
in a single-mode waveguide or flber. Under these circum- 
stances transverse momentum entanglement as well as 
mixed spectral-transverse momentum entanglement are 
automatically suppressed. Thus, for single-transverse- 
mode parametric processes we may limit our analysis of 
the entanglement present to the spectral degree of free- 
dom. It has been shown that a large Schmidt number 
can then be obtained by engineering the photon pair 
generation process to yield the largest possible emission 
bandwidth compatible with the smallest possible pump 
bandwidth [3 . Apart from its role in enhancing the at- 
tainable degree of entanglement, a large generation band- 
width also leads to a small correlation time, deflned as 
the width of the time of emission difference (between sig- 
nal and idler) probability distribution |1H [T2]. A source 
with these characteristics would be useful for applica- 
tions relying on a short time of arrival difference be- 
tween two optical modes, such as quantum optical co- 
herence tomography where the instrument resolution 
is inversely proportional to the correlation time. Let us 
note that an ultra-broadband photon pair source treated 
classically, i.e. where we employ standard non-photon- 
counting detectors, could serve as a substitute for white 
light sources (e.g. based on self-phase modulation in 
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fibers) with one clear advantage: while white light spec- 
tra often have a rather irregular shape, ultra-broadband 
parametric processes can be engineered to have a nearly 
flat spectrum. Likewise, a source with these character- 
istics could lead to parametric amplifiers with an excep- 
tionally large bandwidth. 

The generation of ultra-broadband two-photon states 
is possible with sources based on parametric down- 
conversion (PDC) relying on second-order non-linear 
crystals [131 [121 HSl dZ]- In a recent experiment, it was 
shown that selecting the non-linear crystal and pump 
frequency so that the emitted light is centered at the 
zero group velocity dispersion frequency of the non- 
linear medium leads to PDC light with a remarkably 
broad spectrum; a full width at half maximum gener- 
ation bandwidth of l.OSO^m centered at 1.885/im was 
demonstrated [TB]. Nevertheless, an important limitation 
of an approach based on PDC, is that because broadband 
emission occurs at the zero dispersion frequency, the cen- 
tral emission frequency cannot be freely selected. Note 
that while periodic poling in non- linear crystals [TH] can 
be exploited to yield phase-matching at arbitrary wave- 
lengths, it cannot be used to manipulate the zero disper- 
sion frequency, and hence to select the ultra-broadband 
central frequency; thus, the experiment discussed above 
operates in a spectral region which is inconvenient for 
many applications. This leads to the motivation for the 
present work, in which we aim to develop fiber-based 
ultra-broadband photon pair sources, with far greater 
freedom for selecting the central generation frequency 
and other emission characteristics. 

In this paper we analyze the spectral properties of pho- 
ton pairs generated through the process of spontaneous 
four wave mixing (SFWM) in optical fiber. We concen- 
trate on the case of quasi-monochromatic pumps which 
may be either frequency degenerate or non-degenerate, 
and where all the fields propagate in the fundamental 
fiber mode. While our theory can be applied to any fiber, 
we illustrate our discussion with a numerical exercise for 
the specific case of photonic crystal fiber, composed of 
a solid core and a cladding with an array of air holes. 
The main motivation for using PCF is that the resulting 
dispersion properties can be tailored by careful choice 
of the fiber design parameters including the size, loca- 
tion and shape of the air holes, which in turn permits 
the generation of photon pairs with tailored properties. 
We analyze the conditions which must be imposed on 
the pump frequencies and the fiber parameters to permit 
ultra-broadband photon pair generation. Note that while 
broadband four-wave- mixing in optical fiber has been an- 
alyzed from a classical perspective by several groups f20], 
such an analysis has not been presented for the sponta- 
neous, non-classical case. We also discuss a quantum in- 
terference effect which can occur in the process of SFWM 
when the pump includes multiple spectral lines, for ex- 
ample corresponding to degenerate and non-degenerate 
pumps. 



II. SPONTANEOUS FOUR WAVE MIXING 
THEORY 

We study the spontaneous four wave mixing process in 
single-mode fiber with a third-order nonlinearity x^'^'- I^i 
this process, a photon pair, comprised of one photon in 
the signal mode. Eg, and one photon in the idler mode 
Ei, is created by joint annihilation of two photons from 
the pump fields Ei and £'2- In this paper we focus on 
source geometries which permit the generation of par- 
ticularly broadband signal and idler photon-pair states. 
Remarkably, this is possible even if the pumps are nearly 
monochromatic. Our theory in this paper is valid in the 
limit of monochromatic pumps, and for all generation fre- 
quencies. In contrast, the theory presented in a related 
paper from our group frOj for broadband pumps, is valid 
only for a narrow signal and idler spectral vicinity and 
is therefore unsuitable for the description of broadband 
photon-pair generation. In our present analysis we as- 
sume that all fields are co-polarized and that they prop- 
agate in the fundamental transverse mode of the fiber. 

The quantum state of the generated photon pair in 
an optical fiber of length L can be obtained following a 
standard perturbative approach [21] and is given by 

l^-) = |0),|0)i + K J J duJsduJiF{u;s,u;^)\Lo,)JLu,)^, (1) 

where k is a constant which represents the generation 
efficiency and F^LOs^uJi) is the joint spectral amplitude 
function (JSA), which describes the spectral entangle- 
ment properties of the photon pair 

F{u;s,i^i)^ J duj'ai (w') 02 {uJs + uji - lj') 

X sine [LAk {lo' , ujs, cui) /2] exp [iLAk {uj' , uig, Wi) /2] . 

(2) 

The JSA function is given in terms of the phasemis- 
match function 

Afc {LUl,U!s,Wi) = k (ui) + k {uJs + ~ ^1) ~ k{ujs) 

-A:(c^,)-(7iPi+72P2), (3) 

which includes self / cross-phase modulation contributions 
for the two pumps with peak powers Pi and P2, charac- 
terized by the nonlinear parameters 71 and 72; ai^uj) 
(with i = 1,2) represents the spectral shape of the 
pumps. The energy conservation constraint is appar- 
ent in the argument of the second term of the phase 
mismatch |10j. 

A. Spontaneous four-wave mixing with 
narrowband pumps 

An expression for the JSA (Eq.([2])) in closed analytic 
form, valid for all generation frequencies, can be derived if 
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both pumps are nearly monochromatic. Here we assume 
that the pumps have a Gaussian spectral profile, centered 
at frequencies uji and lu2 , each with a narrow bandwidth 
a. It can be shown that in the limit where both pumps 
are monochromatic, i.e. a —^ 0, the product of the two 
pump envelope functions in Eq.([2]) reduces to 



X 5[uj' - {ids + uJi + wi - LJ2) /2] 



(4) 



The appearance of a Dirac delta function involving lu' 
in Eq.Q allows us to carry out the integral in Eq.(|2]). 
Thus, in the limit of monochromatic pumps, the JSA is 
given by 



AfcW(<5+, S_) = -(71 Pi + 72P2) + Sko + ^ShS+ 



4(2 



Ml 



1 

8(3!) 



Sh6l 



16(4!) 



Sk.Si 



(10) 



where Sk„ 



-r fv_ — 2fc("\ where k'^^ repre- 
sents the nth frequency derivative of k evaluated at 
fl± = ujzd ± (wi — W2)/2, and where fc^"^ represents the 
nth frequency derivative of k evaluated at ojzd- Here, we 
have defined detunings with respect to uJzdi ^^nd further- 
more have defined new variables given by the sum and 
difi'erence of these detunings 



Fcw{^s-, ^i) = N5{uJs +UJi- LUi- LO2) 

X sine [LA kcw {ujs,^i)/ 2] exp [iL A k^w (t^s , ) /2] , (5) 

in terms of a normalization constant N and the phase 
mismatch AA:cm)(ws, w^) which is now a function only of 
Lj^ and LOj 



Akcwiujs,uji) = k [{lOs + UJi +uji- LU2) /2] 
-I- k [{uJs + UJi - LJi+ UJ2) /2] 
-k{ujs)~ k{uji) - (71 Pi + 72P2) . (6) 



Making use of Eqns. [TJ [5][6| we can write down the two 
photon state as 



= 

5+ = Ss ± S,,. 



(11) 



Let us note that the constant term of the expansion, 
(5fco, vanishes if perfect phasematching is achieved for 
pump frequencies uji and uj2 at the generated frequencies 
UJs = uJzd and uJi = uj^d- The JSA can then be written as 



xexp 



iLAfcW(^+,5_)/2 



(12) 



where we have defined 5p+ — 5i + 82- 



\^) = \Q)s\Q)^ + Kjdujf {uj) \uj)s\uj, +UJ2- UJ)^, (7) 

where 



f{uj) = Afsinc[LAfcc„,sm(,(a;)/2] exp[j H/2], 

(8) 

in terms of 



Akcw,sing{^) = Afccu,(cJ,tJi + UJ2 - U)) ^ k{uJi) 
+ k{uJ2) - k{uj) - k{uJl +UJ2-UJ)- (71 Pi + 72P2). (9) 

Note that while \Fcw{uJs,0Ji)\'^ represents the joint spec- 
trum, |/(ci;)p represents the singles spectrum. From 
Eq.(|5]), it is clear that photon-pair generation requires 
the fulfilment of the following two conditions: i) energy 
conservation, or ujs+uji = uji+uj2 and ii) momentum con- 
servation or Akcw ~ 0, with a tolerance which is inversely 
proportional to P. In order to analyze the phase match- 
ing properties near the zero group velocity dispersion fre- 
quency ujzd of the fiber, we express the propagation con- 
stant as a Taylor series centered at uJs = uJi = uj^d- Thus, 
the phasemismatch may be expressed to fourth order as 



B. Step index dispersion model for photonic 
crystal fiber 

Our analysis so far, and the conditions for ultrabroad- 
band photon pair generation to be presented in Sec. |III[ 
are valid for any fiber. In order to carry out specific 
numerical calculations, in this paper we use specialize 
our discussion to photonic crystal fiber. This type of 
fiber consists of a fused silica core surrounded by silica 
cladding with a pattern of air holes which remains con- 
stant along the fiber length. This mixture of air and 
and glass in the cladding results in an average refractive 
index that is considerably lower than that of the core, 
providing a high dielectric contrast, resulting in strong 
optical confinement. This leads to high peak irradiances 
even for modest input powers, which enhances nonlinear 
optical effects such as SFWM. In addition, the dispersion 
characteristics of the PCFs can be engineered by varia- 
tions of the distribution, size and shape of the air holes 
surrounding the core. In particular, it becomes possi- 
ble to choose the zero dispersion frequencies, to tailor 
the SFWM phase-matching properties |22] and to design 
fibers approaching endlessly single-mode behavior |23| . 

The SFWM phase-matching properties are determined 
by the fundamental mode propagation constant, given in 
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terms of the effective refractive index rie// by k{u!) = 
neff{uj)uj/c. We adopt a step-index model, where the 
core has radius r, its index is that of fused silica ns{uj), 
and the cladding index is calculated as ridadi^) — 
f)ns{uj), where / is the air-filling fraction. According to 
Ref. [24j this fiber dispersion model is accurate for 0.1 < 
/ < 0.9. In the context of our work, this model permits a 
straightforward exploration of the spectral entanglement 
properties in {r, /} parameter space. 



C. Phase matching properties for degenerate 
pumps 

Let us now restrict our treatment to the degenerate 
pumps regime, i.e. uji — uj2 = tiJp. In this case, f2+ = 
~ ujzd and therefore, it can be inferred that, except 
for the self/cross-phase modulation term, the series in 
Eq.(10l does not start until third-order terms and all 
terms 6kn vanish. Thus, with Pi = P2 = P and 71 = 
72 = 7, Eq.(10l reduces to 



AfcW(<5+,<5_) = -27P-[M' 



3fc(3) 



4(3!) ^ 8(4!) 



(13) 



For a negligible self /cross-phase modulation contribu- 
tion, an analysis of Eq.(13l reveals the existence of two 
distinct phasematching "branches", where each one cor- 
responds to one of the factors in square brackets vanish- 
ing. We refer to the first as the trivial branch; indeed 
together with energy conservation it may be seen to be 
centered at Ws = oJi = ujp. 

Here we concentrate on the second, non-trivial solu- 
tion. To this end, it is instructive to consider a represen- 
tation of the joint spectrum (see Eq.([5|) in {i5s,(5i} (or 
{ws,Wi}) space. In Fig. [l] we present such a plot for a 
particular fiber (see figure caption), where we have not, in 
contrast with the preceding analysis, resorted to a trun- 
cated Taylor series description. Here, the energy con- 
servation delta function is represented by a straight line 
with negative unit slope, the trivial branch corresponds 
to a line with unit slope and the non-trivial branch is 
represented by a curve which near Sg = 5i — Q is el- 
liptical. The spectral width of each of the two phase 
matching branches is inversely proportional to the fiber 
length (for graphical clarity in our plots we have assumed 
a short fiber) . We have also overlapped a plot of the con- 
tour defined by A/c = 0, shown by the thin blue line, 
which corresponds to frequency pairs which yield perfect 
phase matching. Photon pairs are generated in areas 
of {5s,^i} space where any of the two phase matching 
branches meet the energy conservation locus. 

Fig. Ilia) shows an example of the joint spectral in- 
tensity (cjs, Wi) p for degenerate pumps, obeying ujp > 
Wzd (for which the energy conservation delta function is 
presented by the thick black, solid line), leading to a sin- 
gles spectrum |/(u;)p given by Eq. with three peaks 



as shown in Fig. [T]jb); while the outer peaks correspond 
to the nontrivial branch, the central peak corresponds 
to the trivial branch which is centered at the pump fre- 
quency. Note that if the degenerate pump frequency is 
made equal to oJzd, the two non-trivial peaks merge into 
a single peak, which then overlaps with the trivial peak. 
The long-segment dashed line with negative unit slope, 
represents the resulting energy conservation locus, tan- 
gent to the phasematching curve at ojs = LOi = i^zd- If 
the degenerate pump frequency is tuned down further to 
values such that ojp < oJzd (short-segment dashed line) 
then the non-trivial branch becomes inaccessible, leaving 
only the trivial contribution. Note that if the concav- 
ity of the phasematching ellipse is reversed (see below), 
then spectrally-distinct trivial and non-trivial contribu- 
tions are observed for ujp < cozd instead. 



An analysis of Eq.(13l reveals that the non-trivial 



phase matching branch, to fourth order in the phasemis- 
match, has an elliptical locus in {ujs,uji} space. Note that 
our Taylor series description is valid only in the vicinity 
of the zero group velocity dispersion frequency; for large 
enough detunings from this frequency, the phase match- 
ing locus will deviate from the elliptical shape. Note 
also that for the degenerate pump case, the non-trivial 
phase matching characteristics are independent of the 
pump field. In this case, the phase matching curve is 
tangent at the point given by lOs = LOzd and uoi — tOzd 
(i.e. 6s = 5i — 0) to a line with negative unit slope. 
The curvature at this point is proportional to the ratio 
;j(4)^^(3). jf ratio is positive the concavity is 

oriented towards negative values of (5+, while if the ra- 
tio is negative, the concavity is oriented in the opposite 
direction. 



D. Phase matching properties for non-degenerate 
pumps 

Let us now restrict our attention to a SFWM pro- 
cess in the non-degenerate pumps regime, i.e. uJi =/= lu2. 
For negligible self/ cross phase modulation, an analysis of 
Eq.( 10 ) shows the existence of three distinct phasematch- 



ing branches, two of them trivial and one non-trivial. It is 
straightforward to verify by direct evaluation of Eq. ^ 
in the low pump power limit, that frequency pairs on 



the straight lines 6i — 6s i: 6p- (where 6p- = 61 — 62) 



fulfil perfect phasematching. These are trivial phase- 
matching branches; indeed, together with energy conser- 
vation, they can be shown to be centered at: i)^^ = lui, 
UJi = L02 and \\)lOs = ^^2, — ^i- In order to describe 
the non-trivial branch, we initially restrict our attention 
to a small spectral vicinity around ujzd for which second 
and higher-order terms in Eq. [TO] can be neglected. In 
this case, from Eq. [TO] it may be shown that frequency 
pairs on the line 6i — —6s + 2{2'yP — 6ko)/6ki fulfil per- 
fect phasematching. The fourth order term in Eq. [lO] 
(as well as higher order terms not included in our analy- 
sis) has the effect of adding curvature to this non-trivial 
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FIG. 1: (Color online) Representation of the joint and singles spectrum for a photonic crystal fiber with r — 0.7000/.im and 
/ = 0.9, leading to X^d ~ ^-Kc/uj^d = 0.6335/^m, in the low pump peak power limit and with a fiber length of L = 1cm. Note 
that frequency axes are labeled for convenience with wavelength values, (a) Phasematching function jsinc[LA/Scu)(i^s, '^'01/21^ 
plotted as a function of uJs and Wi; lines with negative unit slope represent energy conservation loci for different values of uip 
around ujzd- Thin blue lines represent the perfect phasematching contour. The joint spectrum is obtained as the intersection 
of the energy conservation locus with the phasematched region in {ujs,LJi} space, (b) Singles spectrum for a pump frequency 
which satisfies uip > uj^d- (c) Similar to (a), in the non-degenerate pumps regime, (d) Singles spectrum for non degenerate 
pumps such that (a;i + 0^2) > Scj^^. 



phasematching branch. 

In order to illustrate this discussion, Fig. ^c) shows 
a representation of the joint spectral intensity, for the 
same fiber as in Fig. [l]ja) , but now in the non-degenerate 
pumps regime (with Ai = 27rc/wi = 0.5400/.tm and 
A2 = 27rc/w2 = 0.7000/im). Note that here we have not 
resorted to a truncated Taylor series description. The en- 
ergy conservation delta function is represented by a thick 
black, straight line with negative unit slope, while the 
trivial branches are represented by two straight lines with 
unit slope. The non-trivial branch is the curved line. As 
in the degenerate pump case, the spectral width of each 
of the three phase matching branches is inversely pro- 
portional to the fiber length. We have also overlapped a 
plot of the contour defined by Afc = 0, shown by the thin 
blue line. Photon pairs are generated in areas of {Ss,Si} 
space where any of the three phase matching branches 
meet the energy conservation locus. 

The joint spectral intensity, represented by Fig. [ijc), 
leads to a singles spectrum with four peaks as shown 



in Fig. [l|d). While the outer peaks correspond to the 
nontrivial branch, the two inner peaks correspond to the 
trivial branches. In the case of non-degenerate pumps, 
depending on the pump configuration, it is possible to 
observe: i)two non-trivial peaks spectrally distinct from 
two trivial peaks (as in Fig. [ijd)), ii)a single non-trivial 
peak spectrally distinct from two trivial peaks (if the 
energy conservation locus meets the non-trivial branch 
tangentially), and iii)two trivial peaks only (if the energy 
conservation locus does not meet the non-trivial branch) . 



E. Effect of self/cross phase modulation 

Let us now consider the effect of the self/cross-phase 
modulation term in the degenerate pumps regime (a sim- 
ilar behavior to that to be described here would be ob- 
served for the non-degenerate pumps regime), which be- 
comes important for a sufficiently high pump peak power. 
For a given pump peak power, this term is an additive 



constant to the phasemismatch at zero power Afcp^Q. 
Thus, a contour diagram in {Ss,6i} space formed by the 
contour Ak = for different power levels, represents a 
contour map of the function Akp^Q. Upon increasing 
the pump peak power, the phase matching contour shifts 
towards areas of {Ss,6i} space characterized by higher 
values of Akp^Q. As an illustration, Fig.[2j^a) shows, for 
a given fiber (see figure caption) , the effect on the phase- 
matching contour of increasing the pump power. In the 
background, dark areas correspond to Afcp=o < and 
light areas correspond to Akp=o > 0. It is clear that as 
the pump power level is increased, the phase matching 
contour shifts from the dark-light interface into the light 
areas. How much it shifts at a particular point on the 
curve depends on how slowly Afcp^g changes with fre- 
quency at that point. In general terms, for frequencies 
generated far from the pump frequency, e.g. on the non- 
trivial branch, this shift is much less pronounced than for 
frequencies generated near the pump frequency. 

For the trivial phasematching branch and degener- 
ate pumps, the self/ cross phase modulation term leads 
to the suppression of perfect phasematching in areas 
of {uJs^uJi} space where Afcp=o < (indicated with 
dark shading in Fig. [2]^ a)). In contrast, in areas char- 
acterized by Akp^o > 0, the self/cross-phase modu- 
lation term splits the trivial branches into two paral- 
lel branches. In Fig. ^h) we show a representation 
of the corresponding joint spectral intensity. The plot 
in shades of gray represents the phasematching func- 
tion |sinc[LAA:cu,(ws, Wi)/2] which is maximized, with 
Akcw — 0, for frequency pairs on the curved black lines; 
the straight black line represents the energy conserva- 
tion locus. Note that while the dark-shaded area in 
Fig. |2|a) can be accessed with a pump in the normal 
dispersion {k'^^^ > 0) regime, the light-shaded area can 
be accessed with a pump in the anomalous dispersion 
(A;*^^^ < 0) regime. The splitting of the trivial branch 
in the anomalous-dispersion regime is a manifestation 
of modulation instability. This phenomenon, which is 
equivalent to four-wave mixing phasematched by the 
self/cross-phase modulation, involves the appearance of 
two spectral sidebands, placed symmetrically around the 
pump [5^. In the time domain, and in the stimulated 
regime, this leads to an ultrafast modulation with period 
An/ 6^1 (where Sil is the spectral separation between the 
two sidebands). 



III. CONDITIONS FOR THE GENERATION OF 
ULTRA-BROADBAND PHOTON PAIRS 

In this section, we derive conditions which guarantee 
the emission of highly broadband photon pairs by sponta- 
neous four wave mixing utilizing narrowband pumps, in 
both the degenerate and non-degenerate pumps regimes. 
We will carry out the derivation of these conditions in the 
low pump power limit; in practice, these conditions are 
applicable for peak pump powers up to a few watts. From 
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FIG. 2: (a) Contours Afcc^ (tJs, 0;^) = for the following peak 
powers: P = 200W, P = 400W, P = 600W, P = 800W and 
P — lOOOW (we have used large peak powers for graphical 
clarity). The background shading represents the sign of the 
phasemismatch in the low pump power limit (Afep=o); the 
dark areas correspond to Afcp^o < 0. (b) Representation of 
the joint spectral intensity for degenerate pumps with a fiber 
length of L = 3cm and a peak pump power of P = lOOOW^; 
we have assumed values for r and / as specified in the caption 
of Fig.[l] 

Eq.(|5| it is clear that SFWM photon pair emission occurs 
at signal and idler frequency pairs {uJs,coi} which satisfy 
energy conservation, i.e. Us + uJi = loi + lu2, and which 
in addition satisfy phasematching, i.e. Akcwi^si^i) — 0- 
In order to obtain a large generation bandwidth, we must 
therefore engineer the phasematched region in the joint 
frequency space {wg, uji} to maximize its overlap with the 
energy conservation locus, which is itself a straight line 
parallel to the S- axis. This translates into a condition 
on the orientation, curvature and location, of the phase- 
matching contour Akcwi^s^i^i) = 0. 

Concretely, we expect broadband photon pair genera- 
tion centered at signal and idler frequencies {usoT^io} if 
four conditions are satisfied by these frequencies: i) they 
satisfy perfect phasematching, i.e. Akcwii^soj^io) = 0, 
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ii) they obey energy conservation, i.e. the energy conser- 
vation locus contains this frequency pair, iii) the phase- 
matching contour, which contains this frequency pair [as 
guaranteed by i)], is oriented parallel to the S- axis, and 
iv) the curvature of the phasematching contour vanishes 
at this frequency pair. Let us note that all degenerate 
frequency pairs ujso = '^io (i-e. which are contained by 
the line 5- = 0) fulfil the orientation requirement iii); 
indeed, from Eq. [lO] it may be shown that the phase- 
matching contour fulfils {d5+ / d5-)\s_=Q = 0. We will 
therefore restrict attention to SFWM geometries centered 
at a degenerate frequency pair. Constrained by 5_ =0, 
the attainment of phasematching [condition i) above] is 
guaranteed if, in addition, (5+ = as can be verified from 
Eq. 10 (note that (J/cq vanishes for 5+ = (5_ = 0). Energy 



conservation [condition ii) above] can then be satisfied for 
(5+ = (5_ = if 5p+ = (see Eqp2|), which in turn leads 
to two possible scenarios involving perfect phasematch- 
ing at LOg — uji — ojzd'- i) degenerate pumps (DP) with 
LOi ~ LJ2 — LOzd (or 5i — 82 = Q), and ii) non-degenerate 
pumps (NDP) with ljJ2 — 2ujzd — ^1, where ^1,1^2 7^ (^zd 
(or Si — —62 with 61,62 0). 

Conditions i) through iii) above can be fulfilled for any 
given fiber; they determine the pump and central genera- 
tion frequencies which lead to the optimum SFWM band- 
width, for a specific fiber. The condition on the curvature 
[iv) above], on the other hand, helps us to select the fiber 
core radius, for a given air-filling fraction (or if our anal- 
ysis were applied to standard step-index fiber, it helps to 
select the core radius for a given index contrast). In order 
to analyze the fulfilment of this condition, it is helpful to 
define C — (1 -l-(5^)^'^/^|, where ' denotes a derivative 
with respect to 6- evaluated at <5_ = 0. This quantity 
represents the curvature of the phasematching contour 
at the origin of the generated frequencies 6-} space. 
It is possible to verify that the curvature associated with 
the non-trivial branch, in the low-power limit, is given by 
C = \k^'^'> /{12k^^'>)\. Thus, the curvature is eliminated at 
radii for which fc'^'*^ ^ k'-^^; in general, as will be dis- 
cussed below, for each /, one value of r exists such that 
= 0. The role of fc'-*) is illustrated in Fig. [s] for the 
degenerate pumps regime, where we present phasematch- 
ing contours, represented by solid lines, in a generated 
frequencies (vertical axis) vs pump frequency (horizontal 
axis) diagram. Here, we express the generated frequen- 
cies as detunings from the pump frequency A = uj — ujp; 
the top half of the diagram corresponds to the signal pho- 
ton, while the bottom half corresponds to the idler pho- 
ton. The three panels show the phasematching contour 
for a PCF with / = 0.5, for three different values of the 
core radius, for which: i)k'-^^ < (with r — 0.9023/im) in 
panel (a), ii)A:'^'*' = (with r = 0.7059/im) in panel (b) 
and iii) /c*-^-* > (with r = 0.5953/im) in panel (c). Note 
that in all three cases, conditions i) through iii) from 
Sec. Ill are fulfilled. It is clear that in panel (b), the 



ating remarkably broadband photon pairs. 

Thus, for a source geometry which satisfies conditions 
i) through iv) above, the phasematching contour is a 
straight line around the zero dispersion point 6^ = 6i = 0, 
which in addition coincides with the energy conserva- 
tion locus. For frequency pairs sufficiently removed from 
6s — 6i — 0, the phasematching contour will depart from 
a straight line due to terms fc'-"-' with n > 5. There- 
fore, the attainable phasematching bandwidth is limited 
by the magnitude of these dispersion coefficients. Let us 
note that a fiber characterized by a lower air filling frac- 
tion / involves a lower core-cladding dielectric contrast, 
and consequently the dispersive effects due to the waveg- 
uide contribution tend to be weakened. Such a weaker 
dispersion leads to reduced dispersion coefficients /c("\ 
which in turn limits the departure from a straight line of 
the phasematching contour. Therefore, lower values of / 
permit greater SFWM bandwidths. 

The process of SFWM leads to a remarkable symme- 
try between the two pumps configurations (DP and NDP) 
discussed above. Note that the two-photon state is deter- 
mined by the phasemismatch function Akcw,singi'^) (see 
Eq. 9]) ; thus, if we can show that the DP and NDP scenar- 
ios lead to identical phasemismatch functions, this would 
imply that the two-photon states in these two scenarios 
are themselves identical. To this end, it is instructive to 
calculate the phasemismatch function difference between 
these two scenarios 



ACsm<,(^) - Afcfj;,,„g(c^) = /C(^i) + k{2Uzd ~ ^1) 

~2k{u;zd) = Ak^J^P^giiVzd-, P = 0) = 6ko. (14) 

As indicated in Eq. |14[ the phasemismatch function 
difference between the DP and NDP cases is equal to 
the NDP phasemismatch, for vanishing pump power and 
evaluated ai uj — Uzdi which is in turn equal to 6kQ (see 
Eq. 10 1. Thus, provided that perfect phasematching is 



phasematching contour has an essentially vertical char- 
acter [due to the fulfilment of condition iv)], indicating 
that even a narrow pump bandwidth is capable of gener- 



achieved at degenerate signal and idler frequencies coin- 
ciding with ujzd, Ak^J^J^„g{uj) = AA:^^^„ (w), and there- 
fore the JSA functions, as given by Eq. [5] for these two 
scenarios are in fact identical. Note that this symmetry 
is lost for sufficiently large pump bandwidths. Note, also, 
that this symmetry applies to the JSA, i.e. to the two- 
photon state, but not to the underlying phasematching 
function. In section [V] we discuss an interference effect, 
which arises due to identical outcomes obtained by two 
indistinguishable pathways, when pumping with DP and 
NDP pumps simultaneously. 

For a given fiber (characterized by the values for r and 
/), the singles spectrum |/(a;)|^ can be computed us- 
ing Eq. |8] for a pump frequency such that conditions i) 
through iii) above are satisfied, in either the DP or NDP 
regimes. Given the symmetry discussed above, these two 
regimes give identical SFWM joint, as well as singles, 
spectra. Solid curves in Fig. Qa) show the resulting op- 
timum full width at half maximum singles bandwidth as a 
function of the core radius, where each curve is computed 
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FIG. 3: (Color online) Phasematching [Akcm{i^s,iOi) — 0] contours for fibers with / = 0.5 and: (a) r — 0.9023pim, (b) 
r = 0.7059/^m and (c) r = 0.5953^m. While the thin red curves represent the phasematching contours for a peak power 
of P = 5W, the overlayed thick black curves were computed in the low pump power limit. Vertical dashed lines mark the 
zero dispersion frequency on the pump frequency axes. Note that for r = 0.7059/im, the contour is essentially vertical at 
Ap — Xzd ~ 0.7292/im, yielding a wide generation bandwidth of signal and idler photon pairs. The areas shaded in black 
represent non-physical regions (where ujs or coi would have to be negative to satisfy energy conservation). 



for a fixed value of the air-filling fraction. Here we have 
assumed a fiber length of L = 25cm, strictly monochro- 
matic pumps and a vanishingly low pump peak power 
(we consider the effect of realistic pump power levels in 
the next section). Fig. ^h) shows the zero dispersion 
frequency, which if conditions i) through iii) are satisfied 
corresponds to the central SFWM frequency, plotted vs 
the core radius for different values of the air filling frac- 
tion. These plots show that while spectra centered in 
the infrared tend to be broader than those centered in 
the visible, there is a considerable fiexibility for choos- 
ing the central emission frequency while maintaining a 
large emission bandwidth (as compared to an approach 
based on PDC with x*-^-* crystals). The spectra from 
which the curves in Fig.|4|^a) were calculated typically in- 
clude a main emission mode involving an essentially fiat 
spectrum and in some cases pairs of peaks around the 
main emission mode. For solid- line curves in Fig. [4][a), 
the valleys between satellite peaks and the main emis- 
sion mode do not drop below 0.5 of the maximum value. 
For certain radii, however, the spectrum exhibits pairs 
of well-defined satellite peaks around the main emission 
mode, resulting in voids in the spectrum. The bandwidth 
indicated in Fig. |4|a) by dashed-line curves is calculated 
from the outer-most slopes associated with these satellite 
peaks. It is important to point out that in the case of the 
NDP regime, these satellite peaks are in fact centered at 
the pump frequencies, and are due to trivial phasematch- 
ing. Fig. [5] shows a curve similar to that in Fig. Qa) for 
/ = 0.1 including the complete branch associated with 
satellite peaks. As can be seen, there are fiber core radii 
for which there are in fact more than one pair of satel- 



lite peaks. This is clear, for example, for the core radius 
labelled as C, for which there are two pairs of satellite 
peaks. 

In Fig. |4|a), each solid- line curve is interrupted to- 
wards smaller radii, at the radius for which the zero dis- 
persion frequency no longer exists [conditions i) through 
iii) assume the existence of a zero dispersion frequency]. 
It is apparent from the figure that the generation band- 
width tends to be larger for smaller air filling fractions. 
It is also apparent that for a given air-filling fraction, the 
attainable bandwidth exhibits a peak near the core ra- 
dius rf.{4.)^Q for which the condition k'^^^ = is fulfilled 
(shown as red dots). Note that the exact radius yielding 
the maximum bandwidth is slightly larger than rj.(4)^Q. 
This is due to the infiuence of higher-order dispersive 
terms. In particular, it turns out that a small degree 
of curvature, with fc*^^^ ^ 0, can enhance the attainable 
bandwidth, if fc^^^ and fc^^^ have opposite signs. These 
parameters do in fact have opposite signs over a small 
range of radii larger than rj,(4)^oi explaining the shift of 
the peak from r — rj.(4)^g. 

An interrupted dashed-line curve indicates that the 
emitted spectrum reaches the edge of the range of valid- 
ity of the Sellmeier expansion used to compute dispersion 
for fused silica. 

It is possible to re-interpret the phasematching con- 
tours in the low pump peak power limit, of Fig. |3] as 
non-degenerate pump frequencies in the vertical axis and 
degenerate signal and idler central frequencies in the hor- 
izontal axis. This is inverted with respect to the original 
interpretation, i.e. degenerate pumps in the horizontal 
axis and non-degenerate signal and idler in the vertical 
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FIG. 4: (Color online) (a) Solid lines: FWHM bandwidth in 
the main emission mode as a function of r for different val- 
ues of /. Dotted lines: FWHM bandwidth calculated from 
the outer slopes of satellite peaks; these lines are shown in- 
terrupted when the emitted spectrum reaches the edge of the 
range of validity of the Sellmeier expansion for fused silica. 
For each /, the red dot corresponds to the radius for which 
^(4) _ Q fulfilled, (b) Zero dispersion frequency vs core 
radius r, for different values of the air filling fraction /. Note 
that when conditions i) through iii) of Sec. Ill are satisfied, the 
zero dispersion frequency coincides with the central SFWM 
emission frequency. 



axis. The NDP pump frequencies which satisfy condi- 
tions i) through iii) above can now be determined by set- 
ting the degenerate signal and idler frequency at the zero 
dispersion frequency (denoted in the figure by a vertical 
dashed line), and reading out, on the vertical axis, the 
required NDP pump frequencies from the intersection of 
this vertical line with the phasematching contour. Thus, 
the example in Fig. [sf^a) clearly shows the DP pump fre- 
quency coinciding with w^jj and also shows one set of 
NDP pump frequencies. Fig. [sj^b) shows that source de- 



signs for which the phasematching contour curvature ap- 
proaches zero (fc*^**) K, 0), which yields a large generation 
bandwidth in the DP regime, permit an infinite number 
of pump frequencies which can function as NDP pump 
frequencies (throughout the portion of the contour which 
is essentially vertical). This leads to the remarkable con- 
clusion that an identical two-photon state is expected 
for any pair of pump frequencies amongst the possibly 
infinite set comprised of: i)the DP pump frequency, ii) 
possibly one or more discrete NDP frequency pairs, and 
iii) possibly a continuum of NDP pump frequency pairs. 

A key consideration in the design of fiber-based two- 
photon sources is possible contamination from sponta- 
neous Raman scattering, which occurs over a bandwidth 
of ~ 40THz to the red from each pump spectral band. 
Our treatment in this paper does not take into consider- 
ation Raman gain; within the Raman bandwidth, a full 
analysis must take into consideration the combined ef- 
fects of SFWM and spontaneous Raman scattering and 
therefore our theory is not complete [27]. However, 
since here we concentrate on broadband SFWM, where 
the generation bandwidth is much larger than the Ra- 
man bandwidth, our theory does adequately describe the 
overall two-photon state structure. The need for Raman 
suppression, in order to ensure high-quality signal-idler 
correlations, suggests a useful application for the symme- 
try between the DP and NDP regimes. For a configura- 
tion with an optimized SFWM bandwidth, we may pump 
at a pair of discrete NDP frequencies, which if sufficiently 
removed from the main emission mode, guarantees that 
this main emission mode is free from spontaneous Raman 
scattering. 

Note that longitudinal fluctuations of the fiber prop- 
erties may have an important effect on the proper- 
ties of emitted light, particulary on the attainable 
bandwidth [251 . In practice, this will set a fabrication 
tolerance on the amplitude and characteristic period of 
the fluctuations. 



IV. ULTRA-BROADBAND TWO-PHOTON 
STATES: SPECIFIC EXPERIMENTAL DESIGNS 

In this section we present specific designs of fiber- 
based, ultra-broadband photon-pair sources. As we have 
noted in the previous section, a small air-filling fraction 
tends to enhance the resulting SFWM generation band- 
width. Thus, for the specific examples to be considered 
here, we assume fibers with / = 0.1 (which according to 
Ref. IH] represents the lower bound for the range of valid- 
ity of the step-index, effective-medium dispersion model 
which we have employed). We have selected three dif- 
ferent source designs, involving radii labelled A, B and 
C in Fig. [5j For source A, the core radius is chosen so 
that the condition h^^^ = is fulfilled (corresponding to a 
vanishing phasematching contour curvature). For source 
B, the core radius is chosen so that we obtain the largest 
possible flat-spectrum generation bandwidth possible for 
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FIG. 5: (Color online) This figure is similar to Fig. [4|a), for 
the case / = 0.1. Here we have used a larger plotting range 
to better represent the branch of the generated bandwidth vs 
radius curve corresponding to the satellite peaks. The spectra 
associated with points A, B and C are discussed in Sec. |IV[ 



/ = 0.1 (as discussed above, this maximum does not oc- 
cur at the radius for which fc'"'^ = due to the effect 
of higher-order terms). For source C, we have selected 
the smallest core radius for which the two satellite peaks 
become well defined, i.e. such that the rate of emission 
reaches zero between the main emission band and the 
satellite peaks. In what follows we present emission spec- 
tra for these three cases. In all cases we will assume a 
fiber length of L = 25cm, a pump peak power of P = 5W, 
a nonlinear coefficient of 7 = 70W~^km~^, and a pump 
bandwidth of cr = 50MHz, unless stated otherwise. 

In Figs. [6] and [7] we present the two-photon state 
and singles spectrum obtained for an air-filling fraction 
/ = 0.1 and for a core radius of r = rj,(4)^Q = 1.8162/im 
(which corresponds to point A in Fig. [s]). A number of 
choices for the pump frequencies which satisfy conditions 
i) through iii) above exist, and which therefore permit 
the largest possible generation bandwidth for the fiber 
in question: i)Ai — \i = Xzd =^ '^'^c/^zd = 1.2076/im 
in the DP regime, ii)Ai = 0.7252/im, A2 = 3.6070/im 
(one individual set of NDP frequencies), and iii) any pair 
of frequencies symmetrically displaced from uj^d within 
the range 1.1606/im < A < 1.2586/xm; these consti- 
tute a continuum of NDP frequency pairs. It should be 
stressed that given the symmetry between the DP and 
NDP regimes, all of the above choices of pump frequen- 
cies will result in a two-photon state with the same joint 
spectrum \Fcn,{iOs,u;i)\'^. 

In Fig. [6ja) we present a plot of the phasematching 
function |sinc(LAA:ctD/2)P, in {6+,S-} space, for the DP 
regime. We have overlapped a plot of the perfect phase- 
matching contour Afccu)(^+, <J-) = (for graphical clar- 
ity we have interrupted this contour towards the narrow. 



outer branches). Note that the trivial phasematching 
branch exhibits a power-induced splitting into two par- 
allel branches. Fig.[6]^b) shows the corresponding plot of 
the phasematching function in {S+, 6-} space, plotted for 
the NDP regime (where we have chosen pump frequen- 
cies corresponding to the discrete NDP frequency pair); 
specifically, we have chosen Ai = 2ttc/llJi = 0.7252/xm 
and A2 = 27Tc/uj2 = 3.6070/im. In both cases we have 
indicated with a yellow vertical line the energy conserva- 
tion locus; ideally, for a large generation bandwidth the 
phasematching contour and the energy conservation lo- 
cus should coincide. Remarkably, while the phasematch- 
ing functions in the DP and NDP regimes are consider- 
ably different (as can be appreciated from the plots), the 
resulting two photon states are identical as demanded by 
the symmetry derived in the last section. Indeed, the 
portion of the phasematching function which coincides 
with the energy conservation locus may be seen to be 
identical in these two cases. Fig. |7]^a) shows the joint 
spectral intensity \Fcwii^s,^^i)\^ , plotted as the contour 
corresponding to 0.5 of the maximum value (note that 
the narrow diagonal width should be interpreted only 
schematically; the lines used are much thicker than the 
actual spectral width). As expected, from the symmetry 
discussed above, the joint spectral intensity is identical 
for the DP and NDP regimes. Note that in Fig.[7]^a) there 
is a main emission mode centered at the zero dispersion 
frequency, and there are two satellite peaks, shown cir- 
cled, involving highly non-degenerate frequency pairs. 

The large degree of spectral entanglement evident in 
Fig. ^a,), could in principle be computed through the 
Schmidt number, K. However, because of the very large 
ratio of the large to small diagonal widths, it is difficult 
for the sampling used to suffice in giving a reliable nu- 
merical estimate of K. Our calculation gives us a lower 
bound K > 1.7 x 10^ (where K = 1 indicates a factorable, 
or un-entangled, state). 

Fig. l7tb)shows the singles spectrum, calculated from 
Eq. |8] (and therefore assuming ideal monochromatic 
pumps), which once again is identical for the DP and 
NDP regimes. This spectrum exhibits a main emission 
mode centered at Xzd = '^T^c/uzd — 1.2076/im, with a re- 
markably large FWHM bandwidth of 712.2nm. The sin- 
gles spectrum also exhibits two satellite peaks centered at 
Ai = 27rc/wi = 0.7252/im and A2 = 27rc/w2 = 3.6070/im, 
which in the NDP regime with two discrete peaks [sce- 
nario ii) above] are centered at the pump wavelengths. 
The emission bandwidth calculated from the outer slopes 
of these satellite peaks is 2884. 8nm. The fractional gen- 
eration bandwidth for the main emission mode Auj/uJc 
where Aw is the FWHM bandwidth and ujc is the central 
frequency gives a value of 0.546. This large generation 
bandwidth leads to a small correlation time, of r = 3.4fs, 
calculated as the width of the signal/idler emission time 
difference distribution. 

In Fig. [8] we present the singles spectrum obtained 
for an air filling fraction of / = 0.1 and a core radius 
of r = 1.8402/im (corresponding to point B in Fig. |5] 



11 




COS (10 

0.581 0.963 1.361 





,15 



rad/s) 

1.765 2.139 2.507 



FIG. 6: (Color online) Representation of the joint spec- 
tral intensity in {5+,S-} space obtained for a fiber corre- 
sponding to point A in Fig. [5] Here we present a plot of 
\amc[LAkcw{i^a,i^i)/2]\^ (in shades of gray), and we also in- 
dicate the energy conservation locus as a vertical yellow line. 
While (a) corresponds to the DP regime, (b) corresponds to 
the NDP regime. 



which yields the maximum possible flat-spectrum band- 
width for this air-filling fraction) . As in the previous ex- 
ample, a number of choices for pump frequencies which 
satisfy conditions i) through iii) above exist. These in- 
clude i) Ai = A2 = Xzd = = 1.1987/im in the 
DP regime, two discrete sets of NDP frequencies, at ii) 
Ai = 0.8683Aim, A2 = 1.9349/xm and iii) Ai = 0.7306/im, 
A2 = 3.3367/im, as well as iv) any pair of frequen- 
cies symmetrically displaced from LOzd within the range 
1.1835/im ^ A < 1.2145/im; these constitute a contin- 
uum of NDP frequency pairs. Pump frequencies for sce- 
narios i) and ii) are indicated in the figure by vertical 
dashed lines. Note that the slight departure from a flat 
central emission-mode is power-induced (an essentially 
flat spectrum would be recovered for a lower pump peak 
power). This source leads to a remarkably large gener- 
ation bandwidth of AA = 1142. 3nm, including a central 
mode and two adjacent satellite peaks; this corresponds 
to a fractional bandwidth of Alu/oJc — 0.794. Defining 
the bandwidth according to the outer slopes of the outer- 
most satellite peaks, the generation bandwidth becomes 
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FIG. 7: (Color online) (a) Representation of the joint spec- 
tral intensity in {ujs,uji} space, for the same parameters as in 
Fig. [6] Note that there are two highly non-degenerate satel- 
lite peaks, shown circled, (b) Singles spectrum corresponding 
to the joint spectral intensity in (a). 



2612. 2nm. Note that if non-degenerate pump frequen- 
cies are made to coincide with the first set of satellite 
peaks [scenario ii) above], then the main emission mode 
remains free from spontaneous Raman scattering; indeed 
the red rectangles in the figure indicate the ^ 40THz 
spontaneous Raman scattering bandwidth. The large 
bandwidth of the main emission mode leads to a short 
correlation time of r = 2.5fs. 

In Fig. [9] we present the singles spectrum obtained 
for an air filling fraction of / = 0.1 and a core radius 
of r = 1.8471/im (corresponding to point C in Fig. |5|, 
which represents the smallest core radius for which the 
two satellite peaks become well defined). As in the pre- 
vious examples, a number of choices for pump frequen- 
cies which satisfy conditions i) through iii) above exist. 
These include i) Ai — A2 = Xzd = '^T^c/^zd = 1.1964/im 
in the DP regime, two discrete sets of NDP frequen- 
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FIG. 8: (Color online) Singles spectrum corresponding to 
point B in Fig. [5] The red rectangles represent the ~ 40THz 
spontaneous Raman scattering bandwidth centered at the in- 
ner NDP frequency pair. 



FIG. 9: (Color online) Singles spectrum corresponding to 
point C in Fig. [5] The red rectangles represent the ~ 40TIIz 
spontaneous Raman scattering bandwidth centered at the in- 
ner NDP frequency pair. 



cies, at ii) Ai = 0.8347/im, A2 = 2.1112/im and iii) 
Ai = 0.7347/^m, A2 — 3.2195/Ltm, as well as iv) any pair 
of frequencies symmetrically displaced from ujzd within 
the range 1.1821/im < A < 1.2112yLtm; these constitute 
a continuum of NDP frequency pairs. Pump frequencies 
for scenarios i) and ii) are indicated in the figure by verti- 
cal dashed hues. As in case B, the slight departure from 
a flat central emission-mode is power-induced, and would 
be eliminated in the low pump power limit. This source 
leads to a generation bandwidth of AA — 630. 6nm in 
the central emission mode, corresponding to a fractional 
bandwidth of Auj/uJc = 0.495 and to a correlation time 
of T = 4.1fs. Defining the bandwidth according to the 
outer slopes of the inner-most satellite peaks, the gener- 
ation bandwidth becomes 1330. 8nm. Likewise, defining 
the bandwidth according to the outer slopes of the outer- 
most satellite peaks, the generation bandwidth becomes 
2493.6nm. 



as a continuum of NDP frequencies centered at the zero 
dispersion frequency. In this section we explore the effect 
of pumping simultaneously at the frequencies correspond- 
ing to two elements of this set, where coherence between 
the various frequencies involved is assumed to exist. In 
the discussion which follows, we assume that these two 
elements correspond to DP and NDP frequencies. Be- 
cause a specific photon pair can then be created by any 
of two indistinguishable pathways, we expect an inter- 
ference effect to occur when we vary the relative phase 
between these two sets of pump frequencies. This in- 
terference is such that the two individual pathways can 
interfere destructively leading to the suppression of pho- 
ton pair emission, or they can interfere constructively. 
Variation of the relative phase between the two sets of 
pump frequencies, then leads to a sinusoidal modulation 
of the resulting rate of emission. 



V. QUANTUM INTERFERENCE IN 
SPONTANEOUS FOUR- WAVE MIXING 



In section [nil we concluded that SFWM source designs 
exist where the pump frequencies may be chosen from a 
certain set, which may be finite or infinite, in such a way 
that any choice within this set yields an identical two- 
photon state. This set of pump frequencies is composed 
of all those which lead to the fulfilment of conditions i) 
through iii) of section III It includes one DP frequency, 
coinciding with the zero dispersion frequency, and can 
also include individual pairs of NDP frequencies, as well 



In order to analyze this interference effect, we assume 
a pump field including three narrow spectral bands (e.g. 
with a rectangular shape) centered at each of the DP and 
NDP frequencies. It is straightforward to verify that the 
total peak power available in each of the DP and NDP 
regimes should be identical, so that each pump configu- 
ration leads not only to the same joint spectrum, but also 
to the same rate of emission. One way of attaining this 
condition is if all pump spectral bands gave the same 
spectral width, where in addition the DP amplitude is 
\/2 higher than those for the NDP pumps. Specifically, 
we assume that the pump spectrum can be written as 
follows 
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a([x)) = ao [rect (cjatup^ 1 ; Suj) + A/2e'^rect(ijj£)p; Suj) 
+ rect(wjv_DP,2; Suj)], 



(15) 



where Slu is the bandwidth (assumed to be small) for 
each of the pump spectral components, 9 is the phase 
introduced between the two pump configurations and 
rect(a;; Sx) is unity for —dx/2 < x < bxj^ and zero other- 
wise. In the monochromatic limit where iJw ^ , each of 
the spectral bands may be represented by a delta func- 
tion and the joint spectral amplitude may be computed 
analytically from Eq. [2] yield me 



F(t^,,w,) = N'e'^ cos 0FDpiujs,uJi), 



(16) 



in terms of F]jp(uJs,LUi) which represents the joint spec- 
tral amplitude for either the DP or NDP regimes [which 
are identical if conditions i) through iii) above are satis- 
fied], and where N' represents a normalization constant. 
In Eq. 16 we have omitted a number of terms which are 
negligible, within the spectral range of interest, for a suffi- 
cient spectral separation between the DP and NDP spec- 
tral components. Note that the non-linear nature of the 
photon pair generation process implies that the period as- 
sociated with the interference curve is tt. Note also that 
for a sufficient departure from monochromatic pumps in a 
realistic experimental situation, the joint spectral ampli- 
tude may be computed by numerical integration of Eq. |2] 
We have found that for pump bandwidths up to ^ 0.7nm, 
Eq. [T6]may be used. 

As an illustration, we will consider a specific exam- 
ple based on a photonic crystal fiber with core radius 
r = 0.8658^m and air filling fraction / = 0.4. This 
fiber has a zero dispersion frequency such that Xzd — 
2TTc/uJzd = 0.8089^m. Conditions i) through iii) from 
section Hn] are fulfilled for the DP pump frequency with 
Xdp — 2'kc/lo£,p — 0.8089/im as well as for one NDP fre- 
quency pair with Xndp.i — 2itc/u!ndp,i — 0.7904/im and 
Xndp,2 = 2nc/uNDP,2 = 0.8283/im. We have designed 
the source so that the three frequencies, which must be 
mutually coherent, could be taken from the broad spec- 
trum of a femtosecond-duration pulse train, e.g. from 
a Ti:sapphire oscillator. Fig. [To[a) shows the singles 
spectrum, computed numerically from Eq. [2] assuming 
a pump spectral amplitude of the form shown in Eq. |15| 
for a number of different phase values 9. Here we have 
assumed a peak pump power of P = 5W and likewise 
we have assumed that the three pump spectral bands 
have the same spectral width (equal to AA = 0.5nm) 
and where the amplitude of the central one is hig her by 
a factor of \/2, compared to the other two. In Fig. 10 'a), 
the DP and two NDP spectral components are indicated 
by vertical dashed lines. It is clear from the figure that 
as is increased from to 7r/2, the height of the spec- 
trum is reduced until emission is suppressed at 9 — Tr/2. 
Fig. [To[b) shows the total emitted flux (given as the in- 
tegrated spectral flux) as a function of 9. This plot ex- 
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FIG. 10: (Color online) (a) Singles spectrum resulting for a 
fiber with r — 0.8658pini and / = 0.4, pumped by three nar- 
row spectral bands centered at the frequencies corresponding 
to the DP and NDP regimes. The curves shown were com- 
puted for different phase difi'erences (shown in black boxes) 
between the two pump configurations, (b) Integrated spec- 
tra, proportional to the total emitted flux, as a function of 
the phase difference between the two pump configurations. 



plicitly shows the expected sinusoidal oscillations in the 
emitted flux. 



VI. CONCLUSIONS 

We have analyzed the spectral properties of photon 
pairs generated by the process of spontaneous four wave 
mixing in single-mode fibers and for narrow-band pumps. 
Exploiting this analysis, we have derived conditions un- 
der which ultra-broadband photon pair generation with 
quasi-monochromatic pumps is possible. The resulting 
two-photon states are highly entangled as quantified by 
the Schmidt number and exhibit a particularly short cor- 
relation time. We have found that, for a given fiber, 
the attainable bandwidth is optimized if perfect phase- 
matching is attained for signal and idler frequencies co- 
inciding with the zero dispersion frequency. It is pos- 
sible to design a photon pair source that satisfies this 
requirement which is based on degenerate, or on non- 
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degenerate pumps. These two regimes, referred to as DP 
and NDP in the paper, lead to identical resulting two- 
photon states, revealing a remarkable symmetry in the 
process of SFWM. We have shown that this symmetry 
leads to a quantum interference effect when the fiber 
is pumped, for example, simultaneously by the pump 
frequencies corresponding to the DP and NDP regimes. 
This symmetry also permits, in the NDP regime, the gen- 
eration of ultra-broadband photon pairs without contam- 
ination due to spontaneous Raman scattering; note that 
this represents a key concern in the design of fiber-based 
photon-pair sources. Although our theory can be applied 
to any fiber, we have focused our discussion on the use of 
photonic crystal fibers, described through a step index, 
effective medium dispersion model. We have shown that 
for a given air filling fraction in the cladding, the SFWM 
bandwidth is optimized for core radii close to that which 
leads to the fulfilment of the condition fc*^^^ = 0. Likewise, 
we have shown numerically that smaller air-filling frac- 



tions lead to greater SFWM band widths; this is due to a 
weakening of waveguide dispersion resulting from a lower 
nucleus-cladding index contrast. We have presented spe- 
cific experimental designs, in some cases leading to over 
lOOOnm of emitted SFWM bandwidth. We expect that 
these results will be useful in the design of photon pair 
sources to be used in the exploration of high-dimensional 
continuous- variable entanglement . 
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